An SO(5) model of p-wave superconductivity and ferromagnetism 
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We study an SO (5) model unifying p-wave superconductivity and ferromagnetism. If only a single 
p-wave pairing wavefunction is involved, the p-wave superconducting (pSC) and ferromagnetic (F) 
states can be unified to form 10- dimensional multiplet of SO (5). The collective modes in pSC and 
F states are studied from this viewpoint in terms of a non-linear sigma model. 
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Proximity of magnetism and superconductivity is one 
of the most important issues in the physics of strongly 
correlated electronic systems. In heavy fermion sys- 
tems HH, organic conductors |J, and high-T c supercon- 
ductors(HTSC) jl, the interplay between antiferromag- 
netism (AF) and anisotropic superconductivities (SC) 
has been studied. Recently, a theory has been proposed 
to unify AF and d-wave SC (dSC) in terms of the SO (5) 
symmetry fa]. In this theoretical framework, the order 
parameters of AF and dSC form a five dimensional super- 
vector ft. The formation of the pseudo-gap is identified 
with the onset of the amplitude of this supervector \n\. 
The SO (5) rotation in this five dimensional space turns 
AF into dSC and vice versa. A prediction of this the- 
ory is a specific collective mode in the superconducting 
state, so-called 7r-mode which corresponds to the fluc- 
tuation towards the AF. The corresponding 7r-operator 
becomes canonically conjugate to the AF spin fluctua- 
tion below the superconducting transition and, hence, 
can be detected by neutron scattering experiments. This 
gives a possible scenario for the 41mey-peak observed in 
YBa 2 Cu06.9 Furthermore, this unification is sup- 

ported by some of the recent numerical studies 
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Let us now turn to p-wave superconductors. Ferromag- 
netic spin fluctuations have long been known to mediate 
the p-wave pairing in 3 He |l5| . A more recent exam- 
ple is Si'2Ru04 pq ] . This compound shows a supercon- 
ducting state below T c ~ 1.5K. Unconventional pair- 
ing is clearly indicated by the high sensitivity of T c to 
nonmagnetic impurities fllSfl , and the absence of Hebel- 
Slichter peak in NQR spectrum |l7j . Theoretically it was 
proposed to be a two-dimensional p-wave superconduc- 
tor |Q , and this was confirmed by 17 O NMR Knight shift 
measurement [ fl9f showing the temperature-independent 
behavior of in-plane susceptibility. Among the five p- 
wave pairing states suggested from group theoretical ar- 
guments ^0|, d(k) = z(k x ±ik y ) is the only order param- 
eter plf compatible with broken time-reversal symme- 
try observed by /iSR experiments p2[ . Interestingly, this 
compound is close to ferromagnetic order in the following 



aspects. First, NMR data indicate the enhancement of 
ferromagnetic spin fluctuations [ ^3| . Second, high pres- 
sure destroys superconductivity and drives the material 
towards ferromagnetism [Q. Above 3 GPa the super- 
conducting state has disappeared and around 8 GPa the 
in-plane resistivity is proportional to T 4 / 3 consistent with 
the behavior dominated by two-dimensional ferromag- 
netic spin fluctuations. Lastly, among related materi- 
als Srn+iRunOan+j, ferromagnetic order is observed for 
n = (2), 3, oo |25|-|27||, and with increasing n the ferro- 
magnetic transition temperature rises. 

Motivated by this, we study in this paper an SO(5) 
model describing the interplay between ferromagnetism 
(F) and p-wave superconductivity (pSC). The charge Q, 
the order parameter of pSC, and the total spin S con- 
stitute a ten-dimensional superspin obeying the closed 
SO (5) algebra. Note that this algebra was found in the 
context of 3 He superfluidity p8|-ptj||. We consider the 
case where the amplitude of this superspin is developed 
and the rotational fluctuations between F and pSC states 
constitute low- lying excitations. 

In the weak-coupling BCS theory one introduces a gap 
function, A Q/3 (k) = - J2k' n .s K*/3 7 <5(k, k')(c_ k ' 7 Ck'i), 
where V(k, k') represents the electron-electron interac- 
tion. In this paper we focus on odd-parity pairing, 
where the gap function can be expressed as A(k) = 
'Yl, ll d ll {k)\a^a y . If we assume that V(k, k') is nonzero 
only when k and k' are close to the Fermi surface, 
V(k, k') depends only on the directions k and k'. Hence 
the order parameter of spin-triplet SC state is 

^(k) = i ^(k,p)( c _ pQCp/5 )(aV^) Q/3 (i) 

where is the odd-parity part of VQt, k'), i.e. 

(k, p) is an odd function in both k and p. The 
summation over p is restricted to the region close to 
the Fermi surface. The lowest order term in the expan- 
sion of V^'(k, p) in terms of spherical harmonics yields 

(k, p) <x k • p , which corresponds to the pSC order 
parameter (/ = 1). 
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We now introduce a new operator 

7i>(k) = 5 , (k,p)c- p , Q! (CT y o- A1 ) a/ 3Cp / 3, 



(2) 



withg(k, p) cx ^*'(k,p). It obviously satisfies (717, (k)) oc 
d^(k). The k-dependence indicates angular dependence 
of the pair wavefunction in the pSC phase. At this point 
we may recall that in SO(5) theory of the HTSC, the pair 
wavefunction is restricted to d x 2 _ y 2 [^J. A similar simpli- 
fication is expected, if the rotational symmetry in the or- 
bital space is reduced to the discrete point group (tetrag- 
onal for Sr2Ru04). Therefore, we consider from now on 
a definite form of the pair wavefunction plf , which can 
be implemented by fixing the k-dependence of V, which 
allows us to suppress the k-dependence as g(p) — g(k, p) 
unless necessary. Therefore, superfluid 3 He is outside the 
scope of this SO(5) theory, since the spatial part of the 
pair wavefunction is not fixed. 

Let us define the SO (5) generators L a b of rotations in 
the F-pSC space as elements of an antisymmetric 5x5- 
matrix 
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\ E P ,a,/3 c l a Kp°vP- Provided |.g(p)| 2 = 1, the genera- 
tors L a i, satisfy SO(5) commutation rules p8[-|30[|: 

[L a b, L cd ] = i(S ac L b d — SadLbc — &b C L a d + ^bdL a c)- (4) 

Considering that g(k, p) should be an odd function of 
k and p, we see that \g\ 2 = 1 exactly holds with the 
choice g(k, p) = sgn(k • p), with fixed value of k. There 
are other choices, e.g. g(p) = (p x ± ip y )/V2, which is 
possibly realized in Sr2Ru04 [|Tjj2^]. 

We now define a superspin n of SO (5), which can be 
conveniently labeled by two indices as n a b- There is some 
freedom to define n a b, and the simplest choice is n a b = 
L a b- Instead, we can also define n a b similarly to L a b, but 
with an extra factor w p p3| inside the summation over p 
in each operator, satisfying w p = w- p . A transformation 
rule of the superspin n a b under SO (5) is determined by 
the commutation relation 

[Lab, n c d] = i(S ac nbd - <5ad«-& c - h c n a d + <Ww)- (5) 

There is an obvious similarity between fl) and (|B|) in- 
dicating that the commutation relation (o) corresponds 
to the ten-dimensional adjoint representation. In con- 
trast, in the SO(5) theory of HTSC §, the superspin is 
five-dimensional, constituting the vector representation. 



^34,^42,-^23 are the spin operators S. Thus, L a b plays 
both the roles of SO(5) generator and superspin. The 
physical interpretation of each element of L a b is the fol- 
lowing: and 7r M are the order parameters of the F and 
pSC phase, respectively, and Q is the electron number. 

Now we construct a non-linear sigma model for the su- 
perspin L assuming that the length \L\ has already been 
developed and only the direction represents low-lying de- 
grees of freedom. This appearance of \L\ occurs, in gen- 
eral, at a temperature higher than T c . At the present 
stage it is not possible to identify clearly the onset of \L\ 
in Sr2Ru04, although certain crossover behaviors have 
been seen in various quantities, e.g. resistance |[(|. Note, 
however, that amplitude fluctuations of the order param- 
eter are irrelevant in the sense of renormalization group 
in two dimensions. The detailed analysis of the onset of 
\L\ is left for a future problem. 

The Hamiltonian consists of the anisotropy energy in 
the superspin space and the elastic energy, i.e. spatial 
rigidity, as 
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where the generalized velocity v a b — X)c(^ ac ^^ bc ~ 
Lbc^Lac)- For simplicity, we set p a b = Pba, Xab = Xba- 
Gauge invariance and spin rotational invariance fix the 
form of the first term in (^), because neither linear terms 
nor cross terms in L a b are allowed. Since L is tcn- 
dimensional, it might be possible to include all the ro- 
tations in the ten-dimensional space into the first term, 
i.e. SO(10) generators. Nevertheless, we consider here 
only SO (5) generators, because it is the minimal symme- 
try containing both SO (3) spin-rotational symmetry of 
F phase and U(l)-gauge symmetry of pSC. We do not 
know so far whether a construction of SO(10) genera- 
tors in a second-quantized form is possible. It is natural 
to assume that Xab (Pab) in @ are identical with each 
other within the same phase. Thus, we set \ia = Xa5 = 
Xtt (a = 2,3,4), xi5 = XQ, X23 = X24 = X34 = Xs and 
similarly for p. We shall also add by hand one more 
symmetry-breaking term 



H„ = g{n x 7T* - 7T* x tt) 2 (g < 0), 



(7) 



Here we shall set uv 



1, i.e. n n 



L„b to assure that 



which is proportional to (d x d*) 2 . Since the d- vector 
d(k) — z(k x ± iky) realized in Sr2Ru04 corresponds to a 
unitary state, we assume g < to assign lower energy to 
unitary states (7? |j 7?*) in pSC than nonunitary states. 
This kind of term has appeared previously in Ref. |3^] 
in connection with the discussion of the stability of the 
ABM state in a certain region of the phase diagram of 
3 He. In contrast to 3 He, where the continuous rotational 
symmetry of pair wavefunctions admits five terms of this 
kind, here the quenched orbital part allows for only two 
terms: and H' = (\Rcd\ 2 + |Imdj 2 ) 2 . The latter term 
contributes to the anisotropy between F and pSC states. 
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However this anisotropy has been already taken into ac- 
count in anisotropies of x's in the first term of eq.(12) 
since we are considering the case with fixed \L\. There- 
fore, H' does not give qualitatively new effects. 

Let us discuss the parameters in (H). An excita- 
tion associated with the term Q 2 /(2xq) corresponds to 
a plasmon mode with high energy. Therefore, we set 
XQ "C XSiX-k- The ground state is either pSC or F, 
depending on the values of xs and Xw H Xir < Xs 
(Xtt > Xs), the system is in the F (pSC) phase. 

As in [||, one can use Hamiltonian (g) (Q) to study 
collective modes. The equation of motion is given by 

Lab = Q yXXqj - Xbd){L a d, L bd } + - y^XPad ~ Pbd){v a d, 
d d 

+ 2 ^2 Pac^{ v ac, L cd L bd } - 2 X/ Pbo^i v bc, L cd L ad } 

+ 2 Pcd^{v c d, L ca L bd — L cb L ad } + i[H v> L a b], 

c,d 

with {A, B} = AB + BA. Using (||), we can discuss 
collective modes in each phase: pSC and F. 

When Xtt > Xs an d g < 0, all the unitary states form 
a set of degenerate ground states. Without loss of gen- 
erality we consider the state with (£41) = —(L14) = A^ 
and other L ab s are vanishing. When we linearize the 
equation of motion (|J) around this ground state, we get 

L'd = (X^ 1 - Xs 1 )^^ + P ff ^V 2 L c4 (c = 2,3) (9) 

La = - Ps A%V 2 L cl (c=2,3) (10) 

(11) 
(12) 

(13) 



Lastly, the mode in ( |ll| ) and (12) can be neglected, be- 
cause they move up to high energies when long range 
Coulomb interaction is considered. 

When the system is in the F phase (x-k < Xs), we 
linearize (^) around the ground state corresponding to 
(L23} = (S z ) = As and obtain 



L C 3 

L43 
L42 
L c d 



= (Xs ~ Xn )A s L c2 
= (X^ 1 ~ Xs^AsLea 
= Ps A%V 2 L42 
= -psA%V 2 L 



P7r A s W 2 L c2 
p n A s V 2 L c3 



(c=l,5) (15) 
(c=l,5) 
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I'b, 



£51 = p Q AlW 2 L 5 4 

L 54 = (Xq 1 ~ X^)L 5 i ~ PvAlV 2 L bl 

L cd = {c,de {2,3,5}) 



At a first glance ( |l3| ) seems to represent Goldstone 
modes, but it is not the case, because L cd (c, d £ {2, 3, 5}) 
is not coupled to L41. On the other hand, the rotations 
in spin space give rise to Goldstone modes, which are 
described by the coupled equations @ and Jic|) . The 
commutation relation [tt\, S^] = isx^ir^ gives canoni- 
cal conjugate relation between spin and 7r's, when pSC 
order parameter in the r.h.s. is replaced by a nonzero 
expectation value, and L 2 4 = —S y and L34 = S x corre- 
spond to rotations of the d-vector in pSC phase, while 
L21 = tt x + ttI and L31 — n v +irl to rotations toward the 
F phase |32j| . They yield two collective modes with the 
frequency 



io(k) = Alk^[p~s~ix 



-1 
s 



x^ 1 



pvAlk 2 



(14) 



One new aspect here compared with previous discussion 
p2| is that the "momenta" 7r correspond to the fluctu- 
ation toward F state and the frequency is determined 
by the anisotropy energy in superspin space, x^ 1 — x^ 1 . 



(16) 
(17) 
(18) 

= Me {1,4,5}) (19) 

and (jlj) lead to a ferromagnetic spin-wave excita- 
tion with dispersion u> = psA s k 2 , while ( |l5| ) and ( |l6| ) 
couple to give two collective modes, which are fluctu- 
ations toward pSC phase. Anisotropy of Xab between 
pSC and F phase produces a mass in these two modes, 
W= A s (x^ 1 -Xs 1 )+P7rA 3 s k 2 . Finally, @ is not related 
with collective modes, because L cd (c,d € {1,4,5}) does 
not couple with L23. 

Let us comment on the case Xtt = Xs, where the sym- 
metry is higher. At this point unitary pSC states and 
F states are degenerate, and the system undergoes a 
first-order phase transition. All the low-energy collective 
modes have the same dispersion u> — p^Ag ^k 2 . 

In real materials, the spin-orbit coupling is not negligi- 
ble and introduces additional anisotropy for the direction 
of spin and d- vector. In Sr2Ru04, this anisotropy may 
fix the direction of d as d || e 2 f2l| , |22|| . We treat this by 
adding to the Hamiltonian extra terms allowed by sym- 
metry. The combination k x ± ik y gives rise to a finite 
value of (L) || e z , and, hence, with the spin-orbit cou- 
pling AS* • L, we expect terms including S z . The gauge 
invariance allows the following extra terms up to the sec- 
ond order in the order parameters. 

H so = -hS z + h'{S z ) 2 - 2K{(Re7r z ) 2 + (Im7r z ) 2 } (20) 

where k is positive since (7?) oc d \\ e z . These terms alter 
the dispersion ( |l4| ) to 



w ± (fc) = A^^ipsAlk 2 + k) (xs 1 - X* 1 + PrrAlk 2 + k) 



±h 



i.e. the spin-orbit coupling gives a finite mass to these 
modes [|20|. The spin susceptibility is then calculated as 

ImXra = I m Xra oc 2J 22 - s^ a (k)) (21) 



a — ± ; 



=±1 



On the other hand, S z is a constant of motion. Thus, 
double peaks would be observable in neutron scattering, 
due to the spin-orbit coupling. The additional terms in 
H so are approximate in the sense that we restrict the 
Hilbert space of d(k) by fixing the ^-dependence. Strictly 
speaking we have to take into account the tensorial struc- 
ture due to the fc-dependence in the quadratic term of the 
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7r-fLelds, which is introduced by the spin-orbit coupling. 
For the above mode they lead to a minor modification 
only, but leave qualitative aspects unchanged. 

The above SO (5) generators and superspins can be cast 

which is relevant 



into a spinor form as in 1 33 



constructing a microscopic theory. 
4-component spinor 



for 

Let us introduce a 



(c P T,Cpi,.9(p)V pT ,.9(p)*cl pi ), (22) 



where the superscript ' denotes transposition. The re- 
dundancy in \& p is expressed as = — g(p)iZ*_ p with 

or 

-1 



R = 



Anticommutation relations are 



We introduce a representation of Clifford algebra 
{T a , r b } = 25 ab ,(a,b = !,•■■, 5); 











• r 



(2,3,4) _ 



a 
'o* 



r 5 _ / CJV 
\ a y o 



in order to write SO (5) generators and superspins as 
^^E*P rQb *P< "^jE^P^'tp, (23) 



where r afc = -i[r a ,r 6 ]. It is noteworthy that the 5- 
dimensional vector | ^ D \E f l,r ai iI/ p , an analog of the 5- 
dimensional superspin |33|| in the HTSC, vanishes iden- 
tically. That is why the superspin in our SO (5) theory 
is not 5-dimensional. In the SO(5) theory of the HTSC, 
the situation is opposite; the 5-dimensional superspin is 
I E P w P* P ra *P+Q in thc notation of ||, while its 10- 
dimensional counterpart | ^ p w p 4 , p r ab 4 , p+ Q vanishes. 

In conclusion, we have constructed an SO(5) model 
unifying p-wave superconductivity and ferromagnetism, 
which may apply to Sr2Ru04. This model describes 
a phase transition between pSC and F by varying the 
strength of the symmetry-breaking term. It also predicts 
spectra of collective modes, among which are collective 
modes corresponding the fluctuation between pSC and 
F phase. The actual phase diagram close to the transi- 
tion point should be governed by these fluctuations and 
will be discussed elsewhere. It would be interesting if the 
phase transition in Sr 2 Ru04 could be observed by chang- 
ing external parameters such as pressure. However, the 
effect of pressure is less clear in the present model, com- 
pared to SO(5) theory of HTSC, where the doping, i.e. 
change of chemical potential, is directly coupled to one of 
SO (5) generators and gives rise to a supcrspin-flop tran- 
sition. 
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